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ABSTRACT 

Taking into account the rotation of mass-accreting white dwarfs (WDs) whose masses exceed the Chan- 
drasekhar mass, we extend our new single degenerate model for the progenitors of Type la supernovae (SNe 
la), accounting for two types of binary systems, those with a main sequence companion and those with a red- 
giant (RG) companion. We present a mass distribution of WDs exploding as SNe la, where the WD mass ranges 
from 1.38 to 2.3 Mq. These progenitor models are assigned to various types of SNe la. A lower mass range 
of WDs (1 .38 Mq < Mwd < 1 -5 Mq), which are supported by rigid rotation, correspond to normal SNe la. A 
variety of spin-down time may lead to a variation of brightness. A higher mass range of WDs (M WD > 1.5 Mq), 
which are supported by differential rotation, correspond to brighter SNe la such as SN 199 IT In this case, a 
variety of the WD mass may lead to a variation of brightness. We also show the evolutionary states of the 
companion stars at SN la explosions and pose constraints on the unseen companions. In the WD+RG systems, 
in particular, most of the RG companions have evolved to helium/carbon-oxygen WDs in the spin-down phase 
before the SN la explosions. In such a case, we do not expect any prominent signature of the companion imme- 
diately before and after the explosion. We also compare our new models with the recent stringent constraints 
on the unseen progenitors of SNe la such as SN 201 lfe. 

Subject headings: binaries: close — stars: winds, outflows — supernovae: individual (SN 201 lfe) — super- 
novae: general 



1. INTRODUCTION 

Type la supernovae (SNe la) play important roles in astro- 
physics as a standard candle for measuring cosmological dis- 
tances and as main production sites of iron group elements. 
It is commonly agreed that the exploding star is a mass- 
accreting carbon-oxygen (C+O) white dwarf (WD). However, 
it is not clarified yet whether the WD accretes H/He-rich mat- 
ter from its binary companion [single degenerate (SD) sce- 
nario], or two C+O WDs merge rdouble degenerate (DP) 
scenario] (e.g JHillebrandt & Niemeyer 120001: iNomoto et al. I 
l2000l) . 

Observations have provided the following constraints on the 
nature of companion stars. Some evidences support the SD 
model, such as th e presence of cir cumstellar matter (CSM) 
dPatat et al. 1120071 LSternberg etaDl2l)Tll iFolev et al. Mim 
and detections of hydrogen in the circumstellar- inter action 
type SNe (Ia/IIn) like SN 200 2ic (lHamuv et al. I l2003t) and 
PTFllkx (Dild av etalH 12012b . On the other hand, there 
has been no direct indication of the presence of compan- 
ions, e.g., (1) the l ack of compan ion stars in the im- 
ages of SN 201 lfe dLi et al. I 1201 ll), s ome SN la rem- 
nants (SNRs) (ISchaefer & Pagnotta 1120 121) SN 1572 (Tycho) 
dKerzendorfetal. 1 120091) and SN 1006 dKerzendorf et al. I 
|2012j), (2) the lack of ultraviolet (UV) excesses of early-time 
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light curves dKasenll2010l). and (3) th e lack of hydrogen fea- 
tures in the spectra ([Leonard 2007). Both (2) and (3) are 
expected from the collision between ejecta and a companion. 

In particular, detailed observations of SN 201 lfe in M101 
requir e stringent constraints on the progenitor, i.e., iLi et al. I 
(1201 ll) excluded the presence of a red-giant (R G), a helium 
star, or a main-sequence (MS) star of > 3 .5 M Q . iBrown et al. I 
d2012l) further excluded a solar mass MS companion from an 
early UV observation with Swift because of no signature of 
shock interaction between ejecta and a companion. 

The tightest constraints come from no signature of shock 
interaction with the companion. However, if the binary sep- 
aration, a, is much larger than the companion radius, R, 
i.e., a S> R, the solid angle subtended by the companion 
would be much smaller, and so would be th e effect o f shock 
inter action. In their spin-d own scenario, Uusthaml d20TTb 
and iDi Stefano et al. I (1201 lh argued that the donor star in 
the SD model might shrink rapidly before the WD explo- 
sion, because it exhausts its hydrogen-rich envelope before 
the SN la explosion during a long spin-down phase of the 
rapidly rotating, super-Chandrasekhar mass WD. In such a 
case, the companion star is much smaller than its Roche lobe, 
which reduces the shock signature. This also explains the 
lack of hydrogen in the spectra of SNe la and possibly th e 
unseen companion in the SNR (IDi Stefano & Kilicl 120121) . 
lHachisu et al. I (120121) presented possible evolutionary routes 
to super-Chandrasekhar mass WDs in the WD+MS channel 
of the SD model. They also compared the spin-down time 
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of WDs with the companion star's MS lifetime and discussed 
their final states at explosions, although their main-focus was 

to explain the observed extremely lumi nous SNe la. 

In this Letter, we apply our method in lHachisu et al. I (120121) 
to WD+RG binaries, and calculate the WD mass distribution 
beyond the Chandrasekhar mass limit for both the WD+MS 
and WD+RG systems. We then estimate the brightness distri- 
bution of SNe la, assuming that the brightness depends on the 
WD mass. We further confirm that, in most of the WD+RG 
systems, the companion has evolved off from a RG to a he- 
lium (or C+O) WD before the SN la explosion and such a 
compact companion does not show any prominent shock sig- 
natures nor indications of hydrogen. In Section|2] we describe 
our basic assumptions and methods. Section [3] presents our 
numerical results. Section|4]discusses various perspectives on 
the progenitors of SNe la. 

2. BINARY EVOLUTION BEYOND THE CHANDRASEKHAR MASS 

Based on the SD model, we followed binary evolutions in 
which a WD accretes hydrogen-rich matter from its com- 
panion. There are two well studied evolutionary paths to 
SNe la, the WD+MS and WD+RG channels. Our basic as- 
sumption s in binary evo l utions ar e essen tia lly the same as 
those in lHachisu et ail dl999allbl I2008allbl l2012h . Mass- 
accreting WDs blow optically thick winds if the mass transfer 
rate exceeds the critical rate, M cr = 6.68 x 10~ 7 (Mwd/M q - 
0.445)M Q yr" 1 (lHachisu et al. ifl^lNomoto et al. Il200l) . 

The WD winds collide with the secondary's surface and 
strip off its surface layer. If the mass-stripping is efficient 
enough, the mass transfer rate is attenuated and the binary 
avoid the formation of a common envelope even for a rather 
massive secondary of ~ 4-6 Mq. Thus, the mass-stripping 
effect widens the donor mass range of SN la progenitors. We 
have i ncorporated th e mass- stripping effect in the same way 
as in Hach isu et al. I (l2008al 120121) . i.e., the mass stripping 
rate M str i p is proportional to the wind mass-loss rate M w j n d 

as M s trip = ciM w jrid. In this study, we assume c\ = 3 as a 
lower representative valu e for the WD+MS system, because 
lHachisu & Kato I (l2003alfbT) found that c\ is between a few to 
10 to reproduce the optical and X-ray light curve behaviors of 
some supersoft X-ray sources. If we adopt a larger value, we 
could have a more massive secondary. For th e WD+RG sys- 
tems, c i is calculated from Equation (21) of Hachisu et al. 
(fl999ah . 

Mass-accreting WDs spin up beca use of angular mo- 
mentum gain from the accreted matter (lLanger et al. 1 120001: 
lUenishi et al7ll2003t iPiersanti et al. 1120031) . If the WD rotates 
rigidly, its mass can only slightly exceed the Chandrasekhar 
mass of no rotation, Mch = 1.46(Ye/0.5) 2 M Q with electron 
mole number Y e . If the WD rotates differe ntially, however , 
its mass can significan tly exceed Mch (e.g., l Hachisu 1986). 
lYoon & Langer I (120041) concluded that the WD increases its 
mass beyond Mch when the accretion rate to the WD is as 
high as Mwd ^ 10~ 7 M Q yr" 1 . They showed that the gradient 
of angular velocity is kept around the critical value for the dy- 
namical shear instability and that this differential rotation law 
is strong en ough to support WDs whose masses significantly 
exceed M ch (lYoon & Langer 120051) . iPiro I (120081) showed that 
highly differential rotation may not be realized due to baro- 
clinic instability. However, his stability condition is not a suf- 
ficient condition but just a necessary co ndition for instability , 
so that his conclusion is premature (see lHachisu et al.1l2.012l) . 

In the present study, we simply assume that mass-accreting 



TABLE 1 

White dwarf mass distribution of Type Ia 
supernovae" 



WD mass 


WD+MS b 


WD+RG C 


ellipticals' 1 


total 


(Mq) 


(%) 


(%) 


(%) 


(%) 


1.38-1.5 


48.2 


41.1 


19.8 


44.8 


1.5-1.6 


13.7 


22.7 


6.7 


17.9 


1.6-1.7 


13.6 


14.0 


2.0 


13.8 


1.7-1.8 


9.6 


10.0 


0.5 


9.8 


1.8-1.9 


8.2 


6.2 


0.0 


7.2 


1.9-2.0 


2.9 


3.7 


0.0 


3.3 


2.0-2.1 


1.8 


1.6 


0.0 


1.7 


2.1-2.2 


1.3 


0.8 


0.0 


1.1 


2.2-2.3 


0.7 


0.0 


0.0 


0.4 



a metallicity Z = 0.02, binary mass ratio distribution f(q) = 1 

b 'Avd+ms = 0.0035 yr~' in our Galaxy 

c *Avd+rg = 0.0031 yr~' in our Galaxy 

d ^2,0 < 1-2 Mq is assumed for elliptical galaxies 



TABLE 2 
WD MASS VERSUS MAXIMUM 
LUMINOSITY DISTRIBUTION 



WD mass 


ratio 


Am 15 (B) 


ratio" 


(Mq) 


(%) 


(mag) 


(%) 


1.38-1.6 


62.7 


1.1-2.1 


67.4 


1.6-1.8 


23.6 


1.0-1.1 


17.3 


1.8-2.0 


10.5 


0.9-1.0 


10.2 


2.0-2.3 


3.2 


0.7-0.9 


5.1 



il taken from Blond in et al. I<20T2T) 



WDs are supported by the above differential rotation law and 
the mass can increase without carbon being fused at the cen- 
ter as long as Mwd > 1 X 10" 7 M Q yr" 1 (=Mb). We assume 
that, when Mwd < M,, hydrogen shell-burning occurs inter- 
mittently on the WD and the resul tant nova outbursts ejec t 
a large part of the envelope mass (lHachisu & Kato I I200D. 
As a result, the net growth rate of the WD mass is signifi- 
cantly reduced (< 10~ 8 M Q yr" 1 ). Then the timescale of an- 
gular momentum deposition to the WD core would become 
much longer than ~ 10 7 yr and comparable to the timescale 
for the Eddington-Sweet meridional circulation (~ 10 s yr) to 
redistribute angular momentum. This redistribution causes a 
contraction of the WD core and its central density increases 
high enough to trigger an SN Ia explosion before the WD 
mass significantly increases. Thus, if the WD mass increases 
beyond 1.38 Mq, we define the WD mass at the SN Ia ex- 
plosion as the WD mass when the mass transfer rate drops to 

Figures Q] and [2] show the parameter regions that produce 
SNe Ia in the logf-M2 (orbital period-companion mass) 
plane for the metallicity of Z = 0.02. We plot results for four 
initial WD masses of M W d,o = 1.1, 1.0, 0.9, and 0.8 M be- 
cause (1) no C+O WDs of Mwd.o S*, 1 ■ 1 Mq are expected for 
normal metallicity (lHachisu et al. 1120 121) and (2) the 0.7 M Q 
region is too small for the WD+MS systems and none for the 
WD+RG systems. The WDs inside these SN Ia regions (la- 
beled "initial") will increase their masses from (a) Mwd = 1.1, 
(b) 1.0, (c) 0.9, and (d) 0.8 M Q to M WD = 1-38, 1.5, 1.6, 
2.2 Mq (0.1 Mq step from outside to inside contours) and 
reach the regions labeled "final." We stop the binary evo- 
lution when the mass transfer rate decreases to Mb. Both 
the WD+MS and WD+RG systems can produce a super- 
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FIG. 1. — Regions producing SNe la with various initial WD masses are plotted in the logP— Mi (orbital period - donor mass) plane both for the WD+MS 
system (left) and the WD+RG system (right). The initial system inside the region encircled by a solid line (labeled "initial") is increasing its WD mass up to the 
mass of Mwd = 1.38, 1.5, 1.6, 2.1, and 2.2 Mq (from outside to inside) and then reaches the regions labeled "final" when the W D stops growing in mass. 
Currently known positions of recurr ent novae and supersoft X-ray so urces are indicated by a star mark (* ) for U Sco (e.g., Hachisu et al. 2000a), a filled square 
for V Sge (Hachisu & Kato 2003b), a filled triangle for T CrB (e.g., Belczvhski & Mikolaiewska 1998), and an open diamond for RS Oph (e.g. JBrandi et al. I 
12003) . 

Chandrasekhar mass WD up to ~ 2 Mq. 

3. DISTRIBUTION OF WD MASSES AT SN la EXPLOSION 



The SN la regions for different initial WD masses, A/wd o = 
0.7, 0.8, 0.9, 1.0, and 1.1 M Q are calculated both for 
the WD+MS and WD+RG systems. We then estimate the 
SN la birth rate in our Galaxy as ^wd+ms = 0.0035 yr" 1 
and fwD+RG = 0.0031 yr" 1 for the constant star formation 
rate. Here we assume t he initial distrib u tion o f binaries 
given by Equation (1) of Eben & Tutukov] (11984 . i.e., v = 
0.2 J f(q)dq J dM/M 2S J dloga yr" 1 and the distribution of 
mass ratio f (q) = 1 . 

We also estimate the delay time distribution (DTD) of 
SNe la as shown in Figure [3^, where the delay time (fdeiay) 
is the elapsed time from binary birth to explosion. The 
spin-down time of WDs is not included. The computa- 



tional method is the same as that in Hachisu et al. (2008b). 
These values are n ormalized to fit the DTD at 11 Gyr 
dMannucci e t al. 2005). Our DTD shows a featureless power 
law (oc ^deiay) from 0.1 to 12 Gyr, which is consistent with 
Totani et al.'s (2008) observation. The presen t results are es- 
sentially the same as our previous results by Hachisu et al. 
(2008b) in which we assume that the WD explodes as an 
SN la at M WD = 1.38 M . In general, the WD+MS sys- 
tems consist of a young population of SNe la corresponding 
to short delay times (fdeiay ^ 1 Gyr) and the WD+RG systems 
an old population of long delay times (fdeiay ^ 1 Gyr). 

We further calculate the distribution of the WD masses at 
SN la explosions as in Table Q] and in Figure [3j)-d. We also 
plot the number ratio of the WDs for the WD+RG systems 
with the initial companion masses of M2.0 < 1.2 Mfp, which 
are expected to occur in elliptical galaxies (see Hachisu et al. 
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FIG. 2. — Continued from Figure[T] 

2010). It is clear that the WD mass distribution in ellipticals 
is confined into a narrower range of 1.38-1 .6 Mq than in late 
type galaxies. 

4. DISCUSSION 
4.1. Spin-down time and final fate ofWDs 

After the mass transfer rate drops to Mb, the WD stops 
growing in mass. There are three characteristic mass ranges 
for the final evolutio n of WDs toward an SN la explosion (see 
lHachisuetal~ll2012l for detail). 

(1) In the extremely massive case, the differentially rotating 
WD explodes as an SN la soon after the WD mass exceeds 
2.4 Mq owing to a secular instability. This is not the present 
case, because it happens for Mwd.o ^1-2 Mq in only low 
metallicity environments. 

(2) For the mid-mass range of M WD = 1 .5-2.4 Mq, the WD 
is differentially rotating and its mass exceeds the maximum 
mass for rigid rotation. As angular momentum in the WD core 



P (day) 

is lost or redistributed toward rigid rotation, the WD core con- 
tracts until its central density and temperature become high 
enough to ignite carbon. Thus the timescale of contraction 
until the SN explosion is ~ 10 s yr due to angular momentum 
transport by the Eddington-Sweet meridional circulation. 

As for the other an gular momentum transport mechanisms, 
lllkov & Soker I d201 lb showed that magneto-dipole radiation 
leads to spin-down in a typical timescale of ~ 10 s — 10 9 yr for 
Mwd si 1 -6 Mq when the magnetic field of the WD is as weak 
as ~ 10 6 G. They also showed that the r-mode instability is 
not significant in spinning-down WDs. Thus, we here assume 
that the Eddington-Sweet meridional circulation is the most 
effective process for spin-down. 

(3) For the lower mass range of Mwd = 1.38-1.5 M Q , the 
WD can be supported by rigid rotation while it exceeds the 
critical mass of non-rotating WDs for carbon ignition. Thus, 
the WD contracts with the spin-down timescale, which is de- 
termined by angular momentum loss from the WD and thus 
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FIG. 3.— (a) Delay time distribution (DTD) for our WD+MS (blue solid), 
WD+RG (red solid) models, and the total (black solid) of them. The ordi- 
nate is the DTD in units of per century and per 10 10 Lk,p>. Open circles with 
an open box are the observational DTD taken from Totani et al. 1 2008) (for 
< 10 Gyr) and Mannucci et al. 1 2005) (at 1 1 Gyr). Each open box indicates 
a 1 cr error of each measurement, (b) Number ratio of the WD masses for the 
WD+MS systems, (c) Same as (b), but for the WD+RG systems. Hatched re- 
gions indicate the contribution of Mi o < 1 .2 Mq , corresponding to elliptical 
galaxies, (d) Same as (b), but for the total of WD+MS and WD+RG systems. 



depends on the strength of the magnetic field of the WD. The 
final fate of the WD depends on the spin-down time as we dis- 
cuss below. It may take more than ~ 10 9 yr for weak magnetic 
fields of - 10 6 G dIlkov &Sokerll20llh . 

If the spin-down time is not much longer than ~ 10 9 yr, the 
compressional heating due to the spin-down would dominate 
the radiative cooling of the WD (see Equation (7) of lNomotol 
1982). Because the spin-down time is not unique, its variation 
causes a variety of thermal state of WD cores when carbon 
ignites at the center. This would lead to a variation of the 
carbon ignition density and thus a variation of 56 Ni mass and 
brightness of SNe la even for the same WD mass. 

If the spin-down time is much longer than ~ 10 9 yr, on the 
other hand, the central density at the ca rbon ignition could 
become high enough to induce collapse dNomoto & K ondo 
119911) . This collapse produces a quite little amount of 56 Ni as 
~ 10~ 3 M Q (IWanajo e t al. 2009), which might correspond to 
a faint transient. 

4.2. Final fate of WD+MS systems 

In most of the WD+MS systems, the companion remains 
to be an MS (central hydrogen burning) star until the "final" 
stage of evolution. As shown in FiguresQ]and[2] these systems 
have an orbital period shorter than ~ 1 day and a companion 
mass smaller than ~ 2 M©. This type of pre-supe rnova bi- 
naries satisfy the constraint on SN 201 lfe posed by iLi et aTl 
d20H,MMS <3.5M©. 

If M WD = 1.5-2.4 Mq and the spin-down time (~ 10 8 yr) 
is shorter than the MS lifetime of the companion, the WD ex- 
plodes before the companion evolves off the main-sequence. 
In most of these cases, the companion's mass further satisfies 
the condition of M M s < 1 M Q posed by iBrown et aTl (12012b 
on the SNR 201 lfe, because the companion's mass further de- 
creases from the "final" mass in Figures[TJand[2]by the amount 
of roughly AM ~ 10" 8 M© yr" 1 x 10 8 yr~ 1 M©, which is 
transferred to the WD and ejected by nova outbursts. On the 
other hand, if Mwd = 1.38-1.5 M© and the spin-down time 
(~ 10 9 yr) is longer than the MS lifetime of the companion, 
the companion becomes a helium WD and CSM has disap- 
peared. Such a case might corresp ond to the case of no CSM 
like SN 201 lfe dPatat et al. ll20Tll) . 

As a lready discussed in our previous paper (Hachis u et al. I 
2008a), the mass-stripping effect produces a large amount of 
CSM, say, a few to several solar masses (« initial mass minus 
final mass, as seen in Figures [TJand|2|. If the spin-down time 
is short enough, or if the WD is forced to be rigidly rotating 
during accretion due to strong magnetic fields, the WD may 
explode in the CSM. Then we may observe the interaction 
between the ejecta and the CSM like in SNe Ia/IIn. 

4.3. Final fate of WD+RG systems 

In the WD+RG systems, after the mass transfer rate drops 
to the companion RG further evolves and finally be- 
comes a helium (or C+O) WD in a timescale of < 10 8 yr. 
This timescale is shorter than the spin-down time in both 
the mid-mass range (1.5 M© < M WD < 2.4 M©: ~ 10 8 yr 
for the Eddington-Sweet circulation) and lower mass range 
(1.38M© <M WD < 1.5 M©: ~ 10 9 yr for the magneto-dipole 
radiation) WDs. The companion RG has already evolved to 
a WD when the primary WD explodes as an SN la. There- 
fore, the immediate progenitor is a wide binary consisting of 
WD+WD for all the cases. These immediate progenitors sat- 
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isfy all the constraints mentioned in Section Q] It might cor- 
respond to SN 201 lfe, which shows no CSM. It also explains 
the lack of hydrogen in the spectra of SNe l a and possibly the 
unsee n companio ns of SN 1572 (Tycho) dKerzendorf et al. I 
l2009h SN 1006 dKerzendorfetal. 1 120121). and SNR 0509- 
67.5 dSchaefer & Pagnottal l2012l) . 

In the above discussion, we assumed that the spin-down 
time is > 10 s — 10 9 yr. However, if the spin-down time is 
short enough, or if the WD is forced to be rigidly rotating 
during accretion, the WD may explode before the companion 
evolves off the red-giant or asymptotic giant branch. Then we 
may observe the inter action between the eje cta and the CSM 
like in Kepler's SNR dChiotellis et al. 1120 121) and in PTF1 lkx 
(iDildav et al. lHoTl) . 

4.4. Variation of Type la supernovae 

In our progenitor models, various types of SNe la can 
be explained as follows: Normal SNe la correspond to the 
lower mass range of WDs, 1.38 Mq < M WD < 1.5 Mq (or 
< 1.6 Mq). The brightness variation can be explained as 
a variety of spin-down time. The brighter group of SNe la 
such as SN 199 IT correspond to the mid-mass range of WDs, 
Myj D > 1.5 M G (or > 1.6 Mq). We think that for these 
brighter SNe la, the brightness variation stems mainly from 
a variation of WD mass. Here we set the border between 
the normal and brighter SNe at about 1.5 Mq. However, it 
depends on the timescale of the Eddington-Sweet circulation 
which may become longer near rigid rotation. Therefore, this 
border could be as massive as Mwd ~ 1.6 Mq. 



Now we explain the luminosity distribution of SNe la with 
our model. In the observation, peak brightnesses of SNe la 
depend monotonically on the Am\s(B), where Am\s{B) is the 
fi-magnitude decay from the maximum in 15 days. Table [2] 
compares the WD mass distribution in our model with t he ob- 
servational Amis(B) distribution (Blon din et al. I f2012). For 
the brighter group of SNe la (M WD > 1.6 M Q ), the distri- 
bution of Amis(B) is in good agreement with the WD mass 
distribution. This may be a support for our theoretical ex- 
pectation that the brightness of SNe la is determined mainly 
by the WD mass because the thermal state is similar among 
various WD cores due to its relatively short spin-down time. 
For the fainter group of SNe la (Mwd ^1-6 Mq), however, 
the brightness of SNe la depends not only on the WD mass 
but also on the spin-down time as mentioned in the previous 
subsection. 

To summarize, we examined the final fate of the two pro- 
genitor models of SNe la, the WD+MS and WD+RG systems. 
A major part of the WD+MS systems reasonably satisfy the 
stringent constraints on SN 201 lfe in M101. Most cases of 
the WD+RG systems satisf y even more stringent constra ints 
on SNR 0509-67.5 posed bv lSchaefer & Pagnottal d20l2l) . 
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